ABSTRACT An analytical solution for the interaction of cylindrical inhomogeneities in a complex potential was obtained using the analytical element method. This method utilizes expansions in the Laurent and Taylor series for the outer and inner parts, respectively, of each cylindrical inhomogeneity. The availability of individual solutions allowed the application of the superposition principle to obtain the solution of the flow problem in terms of discharge potential. This representation, however, required an iterative method for determining the coefficients of the series expansions, which used the minimum value for the absolute differences between successive iterations and a maximum number of iterations as a stop criterion. The implementation of the algorithm was performed in Python language v 3.6.7 and the numerical results after analysis of the convergence are presented in graphs containing the set of inhomogeneities used and their influence on the behavior of flow lines and potentials. These, through the precision and detail obtained, showed the ability of the MEA to characterize flow and identify preferred channels, evidencing its applicability to real situations.

RESUMO O efeito da interação do conjunto de inomogeneidades cilíndricas no potencial complexo em campo uniforme foi analisado por meio da solução proveniente do método de elementos analíticos. Esta, por sua vez, foi obtida por meio das expansões em Série de Laurent e Série de Taylor nas partes exterior e interior, respectivamente, de cada inomogeneidade cilíndrica. O formalismo do método de elementos analíticos, baseado no princípio da superposição, permitiu obter a solução do escoamento em termos de potencial de descarga. A determinação dos coeficientes complexos direcionou para a escolha de um método iterativo, o qual utilizou o valor mínimo para as diferenças absolutas entre iterações sucessivas e um número máximo de iterações como critério de parada. A respectiva implementação de todos os elementos do artigo foi realizada em linguagem Python v 3.6.7. Quanto à disposição das inomogeneidades, espaçamentos regulares e espaçamentos com posições e raios aleatórios são dois elementos considerados; quanto às condutividades quatro cenários são considerados em cada caso de disposição. Em todos os elementos, funções aleatórias foram utilizadas para caracterizar as condutividades do meio e de cada uma das inomogeneidades. Os resultados numéricos permitem avaliar a influência dos diversos cenários de condutividades, bem como a respectiva interação com os poços de bombeamento e injeção no caso de disposição aleatória. Além disso, a precisão e o nível de detalhamento que podem ser alcançados ilustram a capacidade do MEA para a caracterização do escoamento e identificação de canais preferenciais. Por fim, uma análise da convergência do algoritmo iterativo é apresentada para os cenários avaliados.

Resumo A produção de leite é uma das atividades econômicas mais importantes da agropecuária brasileira. Produtores podem usufruir de programas de melhoramento genético que permitem a identificação dos melhores animais e sua seleção para maximizar a produtividade com a melhor relação custo/benefício, além do aconselhamento do descarte de animais menos produtivos. Objetivou-se estimar a herdabilidade e repetibilidade da produção de leite corrigida para 305 dias (PL305) de um rebanho de bovinos da raça Girolando. Foram analisadas 528 lactações de 251 vacas. Para análise foi definida a distribuição uniforme a priori para efeitos sistemáticos. As distribuições de Wishart gaussiana e invertida foram definidas como distribuições a priori para efeitos aleatórios. Os componentes de variância foram estimados utilizando inferência bayesiana pela função MCMCglmm disponível no pacote MCMCglmm do software R. A convergência foi verificada pelo teste de Geweke disponível no software R. Após a obtenção dos componentes de variância foram estimados a herdabilidade e repetibilidade. A herdabilidade observada foi 0,28, o que sugere que a seleção para esta característica resultará em progresso genético eficiente no rebanho. A maior parte da variância fenotípica é devido a variância ambiental, com isso, o fenótipo dos indivíduos não é um bom indicador da variância genética aditiva. A repetibilidade foi de 0,93, indicando que o primeiro desempenho dos animais é considerado um bom indicador do segundo, podendo ser utilizadas em decisões de descarte.

Abstract Milk production is an important economic activity in Brazil. Dairy farmers would benefit from animal breeding programs that aid in identification and selection of animals with the best cost/benefit ratio to maximize productivity, and additionally provide advice on disposal of less productive animals. This study aims to estimate the heritability and repeatability of milk production corrected for 305 days (PL305) in a herd of Girolando cattle. We analyzed 528 lactations in 251 cows. For the analysis, uniform a priori distribution was defined for systematic effects. Gaussian and inverted Wishart distributions were defined as a priori distributions for random effects. The variance components were estimated based on Bayesian inference using the MCMCglmm function available in the MCMCglmm package of the R software. Convergence was verifed with the Geweke test available in the R software. The heritability and repeatability were estimated from the variance component results. Heritability was at 0.28, suggesting that selection for the milk production trait leads to efficient genetic progress in the herd. Phenotypic variance was mainly due to environmental variance; therefore, the phenotype of individuals should not be considered as indicator for additive genetic variance. Repeatability was at 0.93, indicating that the first performance of the animals based on milk production average is a good indicator of the second, and the data could be used for disposal decisions.

Apresenta-se os resultados de uma pesquisa que teve por objetivo exibir um panorama dos cursos de licenciatura em Física do país a partir de uma análise dos Microdados do Censo da Educação Superior do Instituto Nacional de Estudos e Pesquisas Educacionais Anísio Teixeira (INEP). Justifica-se o estudo pelo seu potencial de contribuição ao gerar subsídios teóricos que possam fundamentar o debate público e acadêmico sobre a emergência de ações voltadas para a melhoria do Ensino de Física na Educação Básica. Analisou-se descritivamente diferentes informações de caracterização dos cursos. Como resultados pode-se verificar a existência de 263 graduações em Física (licenciatura), ofertadas por 153 Instituições de Ensino Superior (IES), bem como, percebeu-se, entre outras coisas, que: a distribuição dos cursos não se dá de maneira uniforme pelo Brasil, centrando-se nas regiões Sudeste e Nordeste; 73% das licenciaturas atualmente em atividade foram criadas depois dos anos 2000; identificou-se uma possível associação entre os indicadores de qualidade do INEP e o desempenho dos egressos em exames como o Exame Nacional de Desempenho de Estudantes (ENADE). Sugere-se que em estudos futuros as grades curriculares dos cursos de Física sejam exploradas com maior detalhamento a fim de analisar pontos de convergência entre elas.

The results of a research that aimed to present an overview of the undergraduate courses in Physics in the country are presented from an analysis of the Microdata of the Census of Higher Education of National Institute of Educational Studies and Research Anísio Teixeira (INEP). The study is justified by its potential contribution by generating theoretical subsidies that can support the public and academic debate about the emergence of actions aimed at improving Physics education in Basic Education. Different information on the characterization of the courses was descriptively analyzed. As results, it can be verified the existence of 263 physics degrees (teacher training) offered by 153 Higher Education Institutions (HEIs), as well as, it was perceived, among other things, that: the distribution of courses does not occur uniformly throughout Brazil, focusing on the Southeast and Northeast regions; 73% of the currently active degrees were created after the 2000s; a possible association was identified between quality indicators of the INEP and the performance of graduates in exams such as National Student Performance Exam. It is suggested that in future studies the curriculum of physics courses be explored in greater detail in order to analyze points of convergence between them.

The problem ofelectromagnetic waves propagation inoverhead transmissionlines has apparently not been solved in a sound manner yet. While the problemdoes not have anexact analyticalsolutionwhen considering the presence of the actual surface of theearth, its approximate solutionintroducingtheoreticalsimplificationsis of formidable practical interest. Using quasi-static approximations Carson obtained integral equations to calculate the electromagnetic field due to a horizontal current carrying wire which is above a lossy ground plane. Carsonhimselfproposedthe first solutionto theseexpressions using power series expansions which does not possess uniform convergence andsincethen there have beenefforts toget a bettersolution. In this sense two clear approaches have been essentially followed. The first one consists on modifying the integrand in such a way that an analytic solution can be obtained. The second one is based on using numerical integration schemes.

El problema de la propagación de ondas electromagnéticas en líneas de transmisión aéreas aún no ha sido resuelto de manera definitiva. Si bien el problema no posee una solución analítica exacta cuando se considera la presencia de la superficie real de la tierra, su solución aproximada, introduciendo simplificaciones teóricas es de gran interés práctico. Usando aproximaciones cuasi-estáticas, en 1926 Carson obtuvo ecuaciones integrales para el cálculo del campo electromagnético generado por la corriente de un conductor horizontal sobre un plano de tierra imperfecto. La primera solución la propone el mismo Carson utilizando expansiones en series, las cuales no poseen convergencia uniforme y desde entonces se han hecho esfuerzos por tener una mejor aproximación. Se han seguido dos enfoques claros, el primero consiste en introducir modificaciones en el integrando de manera que sea posible obtener una solución analítica. El segundo, se basa en la utilización de métodos de integración numérica.

Here we study basic approximation properties with rates of our discrete versions of Picard, Gauss-Weierstrass, Poisson-Cauchy singular operators and of two other discrete operators. We prove uniform convergence of these operators to the unit operator. Also all these operators fulfill the global smoothness preservation property. The discussed operators act on the space of uniformly continuous functions over the real line.

Aquí estudiamos las propiedades de aproximación básica con cocientes de nuestras versiones discretas de operadores singulares de Picard, Gauss-Weierstrass, Poisson-Cauchy y de otros dos operadores discretos. Probamos la convergencia uniforme de estos operadores al operador unitario. Además, todos estos operadores satisfacen la propiedad de preservación de suavidad global. Los operadores discutidos actúan sobre el espacio de funciones uniformemente continua sobre la recta real.

We study the sensitivity of convective precipitation to soil moisture in a continental region in Central Amazon based on Atmospheric General Circulation Model (AGCM) simulations and further analyze the physical processes involved, in particular those related to the diurnal cycle. The region was selected based on the analysis of perpetual January simulations with the University of California at Los Angeles AGCM. It shows a uniform and differentiated behavior in the diurnal cycle of convection, which is strongly associated to surface forcing and in opposite phase to large scale ascending motion. Sensitivity to ground wetness was explored based on a control parameter (which varies between 0 and 1) that multiplies evaporation. Results show that the largest sensitivity appears for relatively dry conditions, confirming previous studies. Although evaporation decreases monotonically as the control parameter is reduced, a slight increase in moisture convergence maintains precipitation almost unchanged up to a value of 0.3 of the parameter. The daily maxima in moisture convergence (excluding the extreme case with no evaporation within the region), convective precipitation, and large scale ascending motion at 500 mb are reached in the simulation with a value of 0.4 of the control parameter. Soil and surface temperature increase with decreasing sensitivity parameter, and show retarded daily maxima and larger diurnal amplitude. The same behavior is observed in the sensible heat flux and planetary boundary layer (PBL) height. These effects combine with the variations in PBL water vapor mixing ratio to produce a gradual sensitivity in PBL top relative humidity, in particularly during the morning hours when the humidification observed in the control simulation is gradually reverted, eventually inhibiting early afternoon convection. However, sensitivity of moist convection to ground wetness can not be fully explained without also considering large scale circulation feedbacks. In the more humid cases, deep convection is responsible for mid tropospheric warming. In the drier cases, convective induced warming is insufficient and a differential large scale subsiding motion is induced in the mid troposphere. In the extreme case, with the control parameter equal to zero, the large scale circulation changes drastically with the development of a shallow direct convergent cell that explains the notorious increase in large scale precipitation.

Se estudia la sensibilidad de la precipitación convectiva a variaciones de la humedad del suelo en una región continental del Amazonas en simulaciones con un modelo de circulación general de la atmósfera (MCGA) y se analizan los procesos físicos que establecen dicha sensibilidad, en particular aquellos relacionados con el ciclo diario. Se utilizó el MCGA de la Universidad de California-Los Angeles (UCLA) para simular condiciones de enero perpetuo. Surge del análisis de dichas simulaciones el interés por una región ubicada en el centro del Amazonas, con un comportamiento uniforme -y diferenciado de las regiones circundantes- en el ciclo diario de la convección, fuertemente asociados a procesos de superficie y en fase opuesta al movimiento ascendente de gran escala. Se eligió dicha región para explorar la sensibilidad a la variación de la humedad del suelo, siendo el parámetro de control un factor multiplicativo (entre 0 y 1) de la evaporación. Los resultados muestran que la mayor sensibilidad se tiene para situaciones relativamente secas, confirmando estudios anteriores. Si bien la evaporación es monótona decreciente con la disminución de la humedad del suelo, el leve aumento de la convergencia de humedad mantiene la precipitación en casi los mismos valores hasta niveles de 0.3 del parámetro multiplicativo. Los máximos de convergencia de humedad (excluyendo el caso extremo de evaporación nula), precipitación convectiva y movimiento ascendente en 500 mb se dan para un valor de 0.4 en el parámetro de control. Las temperaturas del suelo y del aire en superficie aumentan con el parámetro de control, con máximos diarios que se retrasan y ciclos diarios de mayor amplitud. Igual comportamiento muestran el calor sensible y la altura de la capa límite planetaria (CLP). Estos efectos, sumados a la variación de la concentración del vapor de agua, se conjugan en un ciclo diario de la humedad relativa al tope de la CLP que varía gradualmente, sobre todo en su comportamiento durante la mañana, donde acaba revertiendo la tendencia creciente, lo cual inhibe el inicio de la actividad convectiva en las primeras horas de la tarde. La sensibilidad de la convección a la humedad del suelo no puede ser comprendida sin considerar también la retroalimentación de la circulación de gran escala. En los casos húmedos, la convección es responsable principal del calentamiento en la troposfera media. En los casos secos el calentamiento convectivo es insuficiente, generándose un movimiento descendente diferencial cerca de la troposfera media. En el caso extremo la circulación cambia drásticamente, instalándose en capas bajas una celda somera con gran convergencia de humedad que explica el incremento en la precipitación de gran escala.

The article deals with three questions. First, if productive specialization on service industries has shown a uniform pattern at regional level. Second, to what extent productive inequalities amongst regions have been reduced or not and which was the role of services. Finally, what has been the contribution of services to regional economic growth. Taking Spain as the case-study, two hypotheses have been specified to answer these questions. The first is that figures should show an inertial trend of the productive specialization, which has produced a low convergence of regional productive structures. The second suggests that productive specialization has influenced regional economic growth and that services have played an important role in it. Several techniques are used to carry out the analysis. Results confirm to a great extent the hypothesis and make possible to set up some aspects not very well known previously.

¿Ha tenido la especialización productiva en servicios un reflejo uniforme a escala regional? ¿Se ha reducido la desigualdad productiva de las regiones y qué papel desempeñaron los servicios? ¿Cuál ha sido la contribución de los servicios al crecimiento económico de las regiones? Para responder a estas preguntas, tomando el caso de España como referencia, se especifican dos hipótesis. La primera es que el período analizado parecería mostrar una notable inercia en la especialización, lo que determina una baja convergencia de las estructuras productivas. La segunda sugiere que la especialización productiva influye en el crecimiento económico regional y que los servicios han tenido una importante función en este sentido. El análisis utiliza diversos indicadores y técnicas y los resultados confirman en gran medida las hipótesis de partida, aportando detalles y aspectos anteriormente no bien conocidos.

We continue our studies in higher order uniform convergence with rates and in Lp convergence with rates. Namely, in this article we establish some Lipschitz type results for the smooth Picard type singular integral operators and for the smooth Gauss-Weierstrass type singular integral operators.

Continuamos nuestros estudios sobre convergencia uniforme de orden superior con radios y sobre convergencia Lp con radios. Concretamente, en este artículo establecemos algunos resultados de tipo Lipschitz para operadores integrales suves del tipo Picard singulares y para operadores integrales singulares de tipo Gauss-Weierstrass.

O equilíbrio de Hardy-Weinberg é um dos principais assuntos estudados pela Genética de populações. Neste contexto, o presente trabalho aborda a análise e a comparação bayesiana de modelos utilizando o coeficiente de desequilíbrio (D A). Para isso, realizou-se um estudo de simulação no qual as seguintes distribuições a priori foram consideradas: Dirichlet (modelo 1); beta - função degrau uniforme (modelo 2); uniforme - função degrau uniforme (modelo 3); e as prioris independentes uniformes (modelo 4). Exemplos de aplicação a dados reais de grupos raciais também são apresentados e discutidos. As amostras das distribuições marginais a posteriori para os parâmetros de interesse foram obtidas mediante o algoritmo Metropolis-Hastings, o qual foi implementado no software livre R. A convergência das cadeias geradas por este algoritmo foi monitorada pelos critérios de Geweke e Gelman & Rubin, os quais estão implementados no pacote BOA do R. Quanto às comparações entre os modelos, efetuadas por meio do fator de Bayes, observa-se que, para os dados simulados, o modelo 4 é o mais indicado para os casos de D A=0,146, D A=0,02 e D A=-0,02 com n=200; o modelo 2 é o mais indicado para D A=-0,02 e n=50 e o modelo 3 é o mais indicado para D A=-0,02 e n=1000. Para os dados reais, em cada caso analisado, nota-se uma grande diferenciação na escolha de modelos, em que apenas o modelo 1 não é recomendado.

One of the main subjects studied by population genetics is the Hardy-Weinberg equilibrium. In this context, this paper addresses the analysis and comparison of bayesian models used in its evaluation by the coefficient of disequilibrium. For this, it was carried out a simulation study in which the following prior distributions were considered: Dirichlet (model 1), beta - uniform step function (model 2), uniform - uniform step function (model 3) and independent uniform priors (model 4). Examples of application to real data for racial groups are presented and discussed. Samples from the marginal posterior distributions for parameters of interest were obtained by Metropolis-Hastings algorithm, which was implemented in the software R. The convergence of the chains generated by this algorithm was monitored by criteria of Geweke and Gelman & Rubin, which are implemented in the BOA package R. Regarding comparisons between models, performed using the Bayes factor, it was observed that model 4 is the most suitable for the cases of D A=0.146, D A=0.02 and D A=-0.02 with n=200, the model 2 is the most suitable for D A=-0.02 with n=50 and the model 3 is the most suitable for D A=-0.02 and n=1000. For real data, in each case examined, there is a large difference in choice of models, where model 1 is the only one not recommended.

One of the inversion schemes most employed in seismic tomography processing is least squares and derived algorithms, using as input data the vector of first arrivals. A division of the whole space between sources and receivers is performed, constructing a pixel model with its elements of the same size. Spatial coverage is defined, then, as the sum of traveled length by all rays through every pixel that conform the medium considered. It is related, therefore, with the source-receiver's distribution and the form of the domain among them. In cross-hole array, rays do not evenly sample the properties of the medium, leading to non-uniform spatial coverage. It is known that this affects the inversion process. The purpose of this paper, then, was to study the problem of spatial coverage uniformity to obtain travel path matrices leading to inversion algorithms with better convergence. The medium was divided in elements of different size but with an even spatial coverage (named as 'ipixels'), and then it was explored how this improved the inversion process. A theoretical model was implemented with added noise to emulate real data; and then the vector of measured times was generated with known velocity distribution. Afterwards an inversion method using minimum length solution was performed to test the two domain divisions. The results showed that the fact of using ipixels not only improved the inversion scheme used in all cases; but in addition allowed to get convergence where it was impossible to do using pixels; particularly through the method considered. This is a direct result of the improvement of condition number of the associated matrices.

Uno de los esquemas de inversion más empleados en el procesamiento de datos de tomografía sísmica es ei de mínimos cuadrados y algoritmos relacionados, utilizando el vector de primeros arribos como datos de entrada. Se Neva a cabo una division del dominio completo entre emisores y receptores, con el objeto de disenar un modelo de píxeles dei mismo tamano. Se define Ia cobertura espacial como Ia suma de los tiempos de viaje de todos los rayos en cada uno de los píxeles que conforman el medio. Por Io tanto este parâmetro está relacionado con Ia distribución emisor-receptor y con la forma dei dominio entre los mismos. En ei dispositivo cross-hole los rayos no muestrean de igual forma ai medio, conduciendo a una cobertura espacial no uniforme. Se sabe que este inconveniente afecta ei proceso de inversion. El propósito de este artículo fué ei de estudiar ei problema de Ia uniformidad de Ia cobertura espacial a fin de lograr matrices de tiempo de viaje que conduzcan a algoritmos de inversion con mejor convergência. El medio se dividió en elementos de diferente tamano pero con cobertura espacial uniforme (denominados 'ipixels'). Se implemento un modelo teórico con ruido a fin de simular datos reales; y el vector de tiempos se calculo con una distribución conocida de velocidades. Luego se probo Ia convergência de Ias dos formas de division del dominio utilizando ei método de solución por mínima longitud del vector de tiempos. Los resultados demostraron que ei hecho de emplear ipixels no solo mejoró la inversion en todos los casos, sino además permitió lograr convergência en casos donde resultó imposible utilizando pixeles; particularmente con el método utilizado. Este es un resultado directo dei aumento de Ia condición de las matrices asociadas.

Este trabalho tem como objetivo realizar uma análise bayesiana de modelos, por meio do fator de Bayes, para o desequilíbrio de Hardy-Weinberg. Pretende-se também testar a metodologia por meio da simulação de dados e aplicá-la a um conjunto de dados reais. Na definição dos modelos, utilizaram-se as prioris Dirichlet (modelo 1), Beta - função degrau Uniforme (modelo 2), Uniforme - função degrau Uniforme (modelo 3) e as prioris independentes Uniformes (modelo 4) relacionadas aos parâmetros coeficiente de endogamia e proporção alélica. Foi implementado um algoritmo no software livre R para realizar a amostragem pelo Metropolis-Hastings das distribuições condicionais a posteriori dos parâmetros dos modelos. A convergência das cadeias foram monitoradas por meio de procedimentos implementados no pacote BOA do software livre R. As comparações entre os modelos indicaram que o mais adequado, ou seja, o que melhor descreve o fenômeno em estudo, é o modelo 1, em comparação aos demais, tanto para os dados simulados, quanto para os dados reais. Em virtude dos resultados apresentados, pode-se atestar que a abordagem Bayesiana apresentou bons resultados, ou seja, por meio das distribuições a posteriori condicionais completas, foram verificadas a confiabilidade e a precisão da metodologia na comparação dos modelos.

The aim of this research is to perform a Bayesian characterization of the Hardy-Weinberg disequilibrium through the Bayes factor. The methodology is tested by using both simulation study and actual data. It was used the following priors for the Bayesian models: Dirichlet (model 1), beta - step uniform function (model 2), uniform - step uniform function (model 3) and independent uniforms for the inbreeding coefficients and allele frequencies (model 4). Metropolis-Hasting algorithms were implemented using the software R to simulate multiple draws from the posterior distribution. Convergence of the Metropolis-Hasting algorithms was assessed by many methods available at R package BOA. Results showed that the model 1 presents the best performance for both simulation study and actual data. The results also showed that the Bayesian approach provides models that are useful for the analysis of the Hardy-Weinberg disequilibrium and inbreeding coefficient.